Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)-dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. solve hydrodynamic equations on an Eulerian grid using Kurganov-Tadmor algorithm combination new adaptive Runge-Kutta method. Our numerical scheme allows us to start soon after nuclear collision, largely avoiding need integrate it separate pre-equilibrium dynamics module. test code's performance by simulating conformal and non-conformal Bjorken flow as well Gubser flow, whose (0+1)-dimensional solutions are precisely known. Finally, we compare hydrodynamics second-order viscous central Pb+Pb find former's longitudinal profile responds more consistently fluid's gradients along spacetime rapidity direction.